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Abstract. A long-lived and multimode quantum memory is a key component 
needed for the development of quantum communication. Here we present temporally 
multiplexed storage of 5 photonic polarization qubits encoded onto weak coherent 
states in a rare-earth-ion doped crystal. Using spin refocusing techniques we can 
preserve the qubits for more than half a millisecond. The temporal multiplexing allows 
us to increase the effective rate of the experiment by a factor of 5, which emphasizes 
the importance of multimode storage for quantum communication. The fidelity upon 
retrieval is higher than the maximum classical fidelity achievable with qubits encoded 
onto single photons and we show that the memory fidelity is mainly limited by the 
memory signal-to-noise ratio. These results show the viability and versatility of long- 
lived, multimode quantum memories based on rare-earth-ion doped crystals. 


1. Introduction 

Distribution of entanglement over long distances is one of the main challenges of quan¬ 
tum information pp. For quantum communication or distributed quantum computing, 
one of the main pillars is the implementation of quantum memories. Indeed, such a 
device could serve as a light-matter interface between computing nodes and channels 
of communication [2], moreover it is an essential building block of quantum repeaters 
lasiE]- These memories need to be efficient, multimode, long-lived and have high 
hdelities illZli. 

For information technology, a multimode capacity of communication channels 
presents the obvious advantage of increasing the density of information processed (com¬ 
puted or transmitted). In quantum information science, the information is often en¬ 
coded, processed and carried by single quanta, for instance single photons. For quantum 
communication single photons present the drawback of being intrinsically sensitive to 
losses, which implies that any operation will be rather inefficient. One can overcome this 
problem by using multimode encoding. Indeed the multimode capacity of the memory 


Multiplexed on-demand storage of polarization qubits in a crystal 


2 


can lead to multiplexing of operations thus effectively increasing the rate of, for in¬ 
stance quantum repeaters ECO] (i.e. the effective success probability of the operation). 
Depending on the system or the type of operation wanted, multiplexing can either be 
temporal [HI [121 113 ] > spatial [HI [151 [T6| or spectral [T7] . 

Ensemble-based quantum memories couple strongly to light due to their intrinsic 
collective enhancement effect. Rare-earth ions doped crystals are promising candidates 
for quantum memories [IHl [S] thanks to their large bandwidth [20] and long coher¬ 
ence times [2I1E21E3] at cryogenic temperature. They present atomic-like properties 
without the need of heavy experimental apparatus for manipulation. Numerous achieve¬ 
ments in the past years have shown the potential of these lanthanides activated solids 
for quantum memories. With the goal of achieving a functioning and scalable quan¬ 
tum memory with these crystals, several milestones have already been reached: high 
efficiency of 69% |24] . cavity enhanced storage for a preprogrammed delay 123 and for 
spin-wave (i.e. on-demand) storage [26], high multimode capacity [l2l[T3], long storage 
time [2I1122]. Recently the longest coherence time ever (6 hours) was measured [23] in 
the same type of crystal we are using here. On-demand spin-wave storage of classical 
time-bin pulses has been achieved in [27112B] and recently time-bin qubits were success¬ 
fully stored on-demand at the single photon level m- Recently we also demonstrated 
on-demand spin-wave storage of a pulse at the single photon level for about 1 ms, using 
spin-echo techniques to extend the storage time [30]. The storage of polarization qubits 
has however been demonstrated for preprogrammed delays only [3TI l32] l33] . 

Here we show the versatility of spin-wave memory at the single-photon-level by 
storing 5 polarization qubits multiplexed in time for more than half a millisecond (50 
times longer than the previously cited work on qubit storage in crystals). Our device 
is based on an atomic frequency comb (AFC) spin-wave memory in an Eu^’'':Y 2 Si 05 
crystal. We use the inherent temporal multiplexing to AFC memories to boost the rate 
of operation of our device by increasing the number of modes stored. We measure the 
fidelity of the retrieved qubits for several mean numbers of input photons and show that 
it is always above the classical bound of 2/3 for true single photons. Furthermore we 
can model the fidelity of our memory based solely on its phase coherence, efficiency and 
unconditional noise floor. This allows us to conclude that its fidelity is mainly limited 
by its signal-to-noise ratio (SNR). 


2. Memory scheme and experimental setup 

2.1. The AFC memory scheme 

The storage protocol we employ is the atomic frequency comb [3l] quantum memory. 
We here give a brief description of the essential features of the protocol, particularly 
the time sequence of the storage process. For a more complete description we refer to 
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Taking advantage of the optical inhomogeneonsly broadened absorption, we tai¬ 
lor an atomic freqnency comb in the absorption prohle throngh a hne selective optical 
pumping. The peaks of the comb are separated by A which is much smaller than the 
input bandwidth. When a photon is absorbed, the single excitation is delocalized over 
all the atoms. This large collective state will quickly lose its collective coherence due to 
the frequency detuning of each atom with respect to the central frequency of the input 
photon. Since all the participating atoms are located in periodic positions in frequency 
space (i.e. mA with m an integer) they will rephase automatically at a time 1/A and 
thereby cause an echo-type reemission. We can intuitively see the reemission process 
as a temporal diffraction (the inverse Fourier transform of the spectral grating i.e. the 
comb) of the input pulse. As such the protocol is a delay line with a preprogrammed 
reemission time. In order to implement on-demand reemission and which can also give a 
longer storage time, we can transfer the optical coherence created by an absorbed pho¬ 
ton into a spin coherence using a strong vr-pulse called control pulse (see hgure 1(b) and 
(c)). This will freeze the phase evolution due to the tailored inhomogeneous frequency 
prohle. When we want to retrieve the stored input, a second control pulse is applied, 
which will restore the optical coherence and the AFC phase evolution will continue up 
to its completion, with the emission of an echo. This scheme is called an AFC spin-wave 
memory. Note that as shown in hgure 1(c) each temporal modes spend the same time 
in the memory (both in the excited state and in the spin state), such that the storage 
efficiency is the same for each mode. 

As for the inhomogeneonsly broadened optical transition, the spin transition is 
broadened (27 kHz in our case). Using spin-echo techniques we can compensate the 
spin dephasing and thus achieve an order of magnitude longer storage time. The RF- 
sequence used in this work is called an XY-4 sequence [35], |36l |37]. This sequence has 
the particularity of being robust to errors in population inversion, independent of the 
phase of the initial spin state. This is relevant in our experiment since the generation 
of the spin state through optical storage has a fluctuating phase. Indeed in a storage 
experiment at the quantum level, it is crucial to achieve very high quality of refocusing 
therefore avoiding unwanted population in the target spin states, which would cause 
fluorescence (i.e. noise) in the output mode (see [ 30 ] )• Furthermore each of the four 
TT-pulses is chirped in amplitude and frequency (adiabatic pulses [HE]) in order to achieve 
efficient population inversion over the entire 27 kHz bandwidth of the spin ensemble. 

The number of temporal modes that can be stored is given by the number of teeth 
in the comb [Mj- The number of teeth can be increased either by increasing the band¬ 
width (for a hxed 1/A storage time) or by increasing the 1/A storage time (for a hxed 
bandwidth). One has to take extra care of the bandwidth of the input modes, which 
cannot exceed the bandwidth of the AFC. For the full spin-wave storage scheme, room 
has to be made for the transfer pulses and of course these control pulses must have the 
required bandwidth [39] . 


Multiplexed on-demand storage of polarization qubits in a crystal 


4 




(c) 



Figure 1. AFC spin-wave storage of polarization qubits, (a) Experimental setup. 
PBS: polarization beam splitter. BS: beam splitter. HWP: half-wave plate. QWP: 
quarter-wave plate. SPAD: single-photon avalanche photodiode. F-P: Fabry-Perot. 
(b) The hyperfine states of the ground and excited states and the two transitions of 
the chosen A-system in Eu^+:Y 2 Si 05 . (c) Timing of the storage sequence. The total 
storage time is 1/A -I-T 5 , where Tg is the spacing between the control pulses. The RF 
spin-echo sequence is inserted in between the optical control pulses. The RF pulses 
are applied using a coil placed around the crystal with a resonant circuit at 34.5 MHz 
(see (a)). 


2.2. Storing polarization qubits in a birefringent and anisotropically absorbing crystal 

The Eu^’'':Y 2 Si 05 crystal is a birefringent crystal which exhibits anisotropic absorption. 
The crystal is cut along the plane (Dl, D2) and the light is propagating along the b axis, 
where (Dl, D2, b) are the principal dielectric axes. There is a coincidence between the 
principal axes of refraction and absorption so the highest absorption is obtained for light 
polarized along Dl, the lowest for light with polarization parallel to D2 [ 30 ] • One can 
thus decompose an arbitrary polarization state into the orthonormal basis formed by Dl 
and D2. Inspired by [SUES] we are using a unique crystal in double pass conhguration 
and a quarter-wave plate in-between the two passes (see figure [^a)). This ensures that 
both components of the polarization will experience the same absorption and global 
refraction. The polarization is preserved up to the swap operation induced by the 
quarter-wave plate in double pass conhguration. We quantify the remaining anisotropy 
by measuring the optical depth and efficiency for different polarization inputs. We found 
that the optical depth varies of ±15% from the mean value which results in a variation 
of ±9% in the efficiency. 
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2.3. Experimental details 

The Crystal and optical setup. We use an isotopically pure 1cm long europium 
doped yttrium orthosilicate crystal, Eu^''':Y 2 Si 05 , with an europium concentration of 
1000 ppm. We work on the ^Fq —)-®Do transition of site 1 at 580.04 nm whose inhomo¬ 
geneous broadening is approximately 1.6 GHz and the overall absorption coefficient is 
a=2.6 cm“^ [H]. The electronic ground and excited states have hyperhne structures 
[121I1S] which are shown in £gure[^b). The inhomogeneous spin linewidth of the 34.54 
MHz spin transition is 27 kHz |30]. The crystal is cooled down below 4 K in a low- 
vibration cryostat (Montana Instruments Cryostation). 

To produce the necessary yellow light, we use a diode laser at 1160 nm, frequency sta¬ 
bilized on a high hnesse optical cavity under vacuum. The laser is frequency doubled 
and amplified to produce more than 1.8 W at 580 nm. The light is then sent through a 
series of acousto-optic modulators (AOMs) in order to prepare three modes: the input 
mode, the control mode and the stabilization mode. The optical setup is shown in £g- 
ure[^a), where both the control and input modes are represented. The input mode is 
an attenuated coherent state, whose average photon number can be varied and whose 
initial polarization is adjusted with the help of a half and quarter waveplate [Prepara¬ 
tion in £gure[^a)). The input mode is then sent through the crystal in the previously 
mentioned double-pass conhguration. After this it is hltered by a Fabry-Perot cavity 
(hnesse 400 and bandwidth 2.5 MHz) whose resonance is stabilized on the input fre¬ 
quency using the stabilization mode (not shown). The output mode is finally analyzed 
with the help of a quarter and half waveplate followed by a polarization beam splitter 
[Analysis in hgure [^a)). Note here that we are using only one port of the PBS for 
detection. The hltering is necessary due to the high power of the control mode (~ 10^^ 
photons) which creates huorescence and off-resonant coherent emission (see [28]). To 
further enhance the hltering, the input and control modes are not applied in the same 
spatial mode. We detect the output using a single photon avalanche diode (SPAD) 
gated with another AOM (not shown in hgure [^a)) before coupling into a hber with a 
9/rm core. The total transmission from the input of the cryostat to the detector is 7%. 
We measure an SPAD efficiency of 57±3% at 580 nm, with a dark count rate of 15 Hz. 
By measuring the counts at the detector and compensating for the rate of experiment, 
total transmission and detector efficiency we can infer the number of photons at the 
input of the memory. In the following we will always refer to this value. 

AFC preparation. The atomic level scheme is shown in hgure [^b). Through optical 
pumping, we hrst prepare a selected sub-ensemble of ions into the ground state | ± |)g 
[26] . Second we create a 2 MHz wide AFC structure on the | ± |)g I |)e transition 
using a 0.5 s long sequence of precise spectral holeburning [39|. During the whole step 
we optically pump the | ± |)g —)■ | ± |)e transition in order to keep the | ± |)g level 
empty, readying it for spin-wave storage. The entire preparation sequence is described 
in more detail in |2ni EH]- 

Storage sequence The full storage sequence is shown in hgure [^c): hve input modes 
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of 1.25 /is each are sent in the crystal and stored with the help of two control pulses 
of 5 /is duration each. In between the optical control pulses, four population-inverting 
(i.e. TT rotation) RF pulses resonant with the spin transition are applied with a periodic 
spacing of T 5 / 4 . Each pulse has a duration of 120 /rs and has a total frequency chirp 
range of 45 kHz around 34.54 MHz. They are generated with an arbitrary function gen¬ 
erator and sent to a 100 W ampliher. The amplihed signal is then sent through a RF 
circulator, an impedance matched circuit made of two tunable capacitors (see in hgure 
[^a)) and finally a 7-turn coil wrapped around the crystal. The impedance-matching 
circuit is tuned to resonate at 34.54 MHz, allowing us to achieve a Rabi frequency of 47 
kHz. 

We measured a transfer efficiency of 70±5% per optical control pulse. The absorption 
probability of the input mode is measured to be 70±3% (see table in appendix C) for 
the comb setting 1/A = 15 /is, such that the total conversion efficiency from the optical 
mode to the spin-wave is about 50%. The total storage time is l/A-\-Ts = 515 /is. 

In the experiment we repeat this storage sequence Wep = 18 times, with an experiment 
rate of 0.8 Hz, which leads to an average rate of storage sequence of 14 Hz. By proceed¬ 
ing so, we considerably increase the effective rate of the experiment. This repetition 
increases the average population in | ± |)p, which in turn would increase the uncondi¬ 
tional noise probability Pn |2H| with increasing number of repetitions N^ep- To reduce 
this accumulation of population we apply a repump pulse on the | ± |)p-| ± |)e transition 
after each storage cycle. 

In the present experiment, we store and retrieve five temporally multiplexed photonic 
qubits encoded with the same state for more than half a millisecond. The different 
stored states are |H), |V), |D) = ;^(|H)-|-|V)) and |R) = ;^(|H)-|-i|V)) with different av¬ 
erage number of input photons (/i = 0.8, 1.4, 3.6, 8.2). 


3. Theoretical model and experimental results 


We first characterize the fidelity of the memory with bright coherent pulses (/i >> 1). 
This allows us to work in a regime where the noise created by the storage process is 
negligible in comparison with the signal. We can then evaluate the phase coherence of 
the memory by calculating its classical fidelity, in order to determine its effect on the 
total conditional fidelity. This classical fidelity can be inferred from the visibility of a 
classical echo signal: 


W = 




Srn.ir, 


Smax + Sr) 


( 1 ) 


where Smin and Smax are the classical signals measured respectively after an orthogonal 
and parallel analyzer with respect to the input. The fidelity is then given by 


Fc = 


l + w 


Sr) 


Smax + Sj) 


( 2 ) 


2 
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Using pulses prepared in various states and equation we calculate a mean classical 
fidelity of 99.1 ± 0.4 %. This shows the good coherence preservation of the memory. 



Figure 2. Photon counting histogram measured for a |D) encoded input with fi=lA. 
The 5 input and output pulses are shown in blue, starting at 0 and 515 gs, respectively. 
The output is measured in the |D) basis. The noise produced by the memory (the 
unconditional noise floor) is shown in green, which was measured with no input pulse 
but still implementing the entire memory protocol. The CPI and CP2 regions show 
where the 5 gs long control pulses are applied. Note the leakage of the control pulse 
through the AOM gate in region CP2. 

Now that we have shown that the memory preserves the states encoded on bright 
pulses, we lowered the number of photons in the input to investigate the behaviour of 
the memory in the single photon regime. Indeed in this regime the noise produced by 
the memory reduces the visibility of the retrieved state. Figure shows the photon 
counting histogram obtained by storing |D) = ;^(|H) + |V)) with /i=1.4 and measured in 
the |D) basis. Similar temporal histograms are obtained when measuring in the |A), 
|H) (or |V)) and |R) (or |L)) bases. We then calculate the fidelity via tomographic 
state-reconstruction using a maximum likelihood method. The experiment is repeated 
with different mean number of photons /i in the input. In order to assess the quality of 
the multimode capacity of the AFC, we can check the individual conditional fidelities of 
the five stored photonic qubits. The fidelities vary from 83.3 ± 3.8% to 86.6 ± 2.9 % for 
fi=lA, showing that the storage fidelity over the different temporal modes is preserved 
(see table 1^ in appendic A). 

We now investigate how the fidelity varies as a function of the mean photon number 
in the input. We also use the multimode capacity in order to increase the rate of the 
experiment, as it would be the case in a quantum repeater [8]. The counts of the 5 
modes are summed to achieve statistics that would otherwise take 5 times (the number 
of modes stored) longer to accumulate. The results are summarized in table 
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h(%) 

Pn (10 

■^) 

hi 

Fidelity (%) 

0.8 

± 

0.1 

4.3 

± 

0.4 

11.0 

± 

1.0 

0.25 

± 

0.04 

79.5 

± 

0.2 

1.4 

± 

0.1 

3.6 

± 

0.3 

10.1 

± 

1.2 

0.28 

± 

0.04 

85.5 

± 

0.1 

3.6 

± 

0.3 

3.8 

± 

0.2 

10.9 

± 

1.4 

0.29 

± 

0.04 

93.6 

± 

0.1 

8.2 

± 

0.6 

3.7 

± 

0.2 

12.1 

± 

1.4 

0.33 

± 

0.05 

95.7 

± 

0.04 


Table 1. Conditional fidelities measured for various mean input photon numbers p, 
prepared in |D). The fidelity was found via tomographic state-reconstruction using a 
maximum likelihood method, and the errors are estimated via Monte-Carlo simulation. 
Shown are also the storage efficiency y, the unconditional noise floor and the pi 
parameter (see text). 


The measured fidelity is reduced as the mean photon number is lowered. To 
understand this behaviour we derive an equivalent of equation ([^ for the fidelity Fg 
in the case where the stored state is at the single photon level. An additional parameter 
has to be considered in this regime: the unconditional noise floor pn |2B]- It corresponds 
to the probability of recovering a photon in the output time window after storage without 
any input, and allows for the definition of a new parameter pi = Pn/f] where rj is the 
memory efficiency. The pi parameter is the mean photon number in the input that 
gives a SNR of 1 in the ouput. Using this figure of merit, we can estimate the fidelity 
of the memory assuming that the noise of the memory is state independent |29l [30] . As 
shown in table in appendix B this asumption is valid in our case. The fidelity at the 
single photon level can be calculated using equation ([^ 


SI 

rp _ _ max _ 

~ I CQ ’ 

^max ~T ^mir) 


( 3 ) 


where Sf^^^=HT]Fc-\-pn and S'U,^=pp(l—Fc)+Pn, given that the input contains pi photons. 
Then 


r]p,Fc + Pn _ y 

T]pi + 2pn 1 + 2 ^' 


( 4 ) 


This formula is valid for a weak coherent state with a mean photon number p as well as 
for a single photon where /i (with 0 < /i < 1) would be the probability to find the pho¬ 
ton before the memory. In figure we compare the fidelities predicted by equation (|^ 
and the measured ones. The agreement is excellent for all mean photon numbers, and 
underscores that the fidelity is limited by noise produced by the memory. We note that 
a higher efficiency would lead to much higher fidelities for mean photon numbers below 1. 


The fidelities in table are always above the classical limit of 2/3 valid for single¬ 
photon qubit storage H, showing that our device is capable of faithfully storing single 
photons with high SNR. However as it has been shown in [32] and [13], one can also 
characterize the quantum nature of a device using weak coherent states. To that end 
we have to consider the maximum storage fidelity achievable for coherent states with a 
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Figure 3. (a) Conditional fidelity as a function of the mean number of photons in the 
input. The dashed red line represents the maximum achievable fidelity with a classical 
memory and single photons. The black curve (grey area) is a plot of equation Q for 
/ri=0.29 (± 0.04) and Fc=99.1 %. The black dashed curve is the maximum achievable 
fidelity with a classical memory and taking advantage of the non-unit efficiency and the 
non-definite (Poissonian) number of photons intrinsic of coherent states (see appendix 
D for more details), (b) Real part of the process matrix as obtained from the quantum 
process tomography. 


classical storage device. For instance a classical memory could measure the input, try to 
deduce its state based on the so gained information, store the result and hnally prepare 
and resend a new state (a measure-and-prepare strategy). The hdelity it could reach 
for a single-photon qubit is limited to 2/3 [H], but in the case of a coherent state, the 
mimicking device could extract more information for states containing more photons. 
Furthermore such a classical storage device could take advantage from the non-unit 
efficiency of the memory by sending a signal only when the number of photons in the 
input is above a certain threshold, which would enable it to reach even higher hdelities 
[291 03] (see hgure in appendix D). 

Here we additionally consider the fact that our device does not fully absorb the 
input state, thereby producing a transmitted signal in addition to the output signal. 
This means that a classical storage device should in principle reproduce the transmitted 
signal as well. As explained in further detail in appendix D the hdelities attainable 
by the classical counterpart are lower than the measure-and-prepare strategy on the 
output alone. This is because the transmission efficiency is higher than the output 
efficiency. This new criterion is plotted in hgure |^a) (dashed line) and with respect 
to it, our memory can be called quantum for the inputs with ju= 1.4, 3.6, 8.2 but we 
cannot conclude with certainty using the state containing p = 0.8 photons on average in 
the current state of the experiment. Nevertheless we would like to emphasize that the 
memory was developed with the goal of storing single photons, which we proved to be 
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feasible provided that the probability of finding the photon before the memory is higher 
than /ii = 0.29. 

Finally, we fixed the mean nnmber of photons to be p=1.4 and we probed onr 
memory with the states |H), |V), |D) and |R) in order to check if the memory can store 
any of those states in the single photon regime. Two methods were nsed to characterize 
the memory. The first method simply consisted in calcnlating the average fidelity, which 
was found to be 84.1 ± 0.1 %, with low deviation from one state to the other (see table 
in appendix B). The second method consisted in constructing the process matrix y 
associated with our device [16]. It is defined as follows: if we input an arbitrary state 
Pin the memory will output the state pout defined by 

3 

Pout ^ ^ X,kMkPin^l y ( 5 ) 

k,l=0 

where ak,i are the Pauli matrices. 

The memory is perfect if and only if yoo=l and all the others contributions are zero: 
in this case Pout=pin- The real part of the reconstructed y matrix is shown in figure 
[^b) indicating the good robustness of our memory (yoo=0.762) which is limited by the 
SNR. However, phase and bit-flip errors are present in the process {X 22 and X 33 errors), 
which possibly could be explained by a residual anisotropic absorption as a function of 
polarization. Indeed, our method for compensating the anisotropy is not perfect. A 
weak residual anisotropy will render the efficiency of the memory weakly polarization 
dependent hence the SNR will slightly vary from one state to the other. We think that 
with further work this polarization absorption dependence can be reduced. Another 
approach would be to increase the efficiency of the memory hence its fidelity, as already 
mentioned above. 

4. Conclusions and Outlook 

We have shown that the AFC spin-wave memory protocol can store multiple qubits at 
the single photon level. We have shown faithful storage of five temporally multiplexed 
polarization qubits for more than half a millisecond in a rare-earth-ion doped crystal. 
Our work shows that the spin-wave manipulation used for storage time extension does 
not degrade the fidelity of the stored states, paving the way for very long storage times 
of quantum states beyond the spin coherence time through dynamical decoupling. The 
excellent agreement of the measured fidelities with those predicted by our simple model 
allows us to foresee the possibility of very high fidelities (a crucial prerequisite for 
quantum repeaters) if we increase enough the efficiency using for instance an impedance 
matched cavity around our crystal [2S]- 
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APPENDIX A: Results of the individual temporal modes 

Here we take advantage of the multimode capacity of our memory by storing 5 temporal 
modes all encoded with the same qubit (|D)). The following table shows that the 
conditional hdelity is the same for all hve temporal modes: 


mode 


vi%) 

Pn (10 



Fidelity (%) 

1 

1.2 

± 

0.1 

3.5 

± 

0.6 

8.8 

± 

1.3 

0.25 

± 

0.08 

84.9 

± 

3.6 

2 

1.5 

± 

0.1 

4.3 

± 

0.6 

12.0 

± 

1.5 

0.28 

± 

0.07 

86.6 

± 

2.9 

3 

1.5 

± 

0.1 

3.2 

± 

0.5 

9.0 

± 

1.4 

0.28 

± 

0.08 

86.4 

± 

3.5 

4 

1.5 

± 

0.1 

3.5 

± 

0.6 

10.5 

± 

1.4 

0.30 

± 

0.08 

85.7 

± 

3.2 

5 

1.5 

± 

0.1 

2.6 

± 

0.5 

9.4 

± 

1.2 

0.36 

± 

0.10 

83.3 

± 

3.8 


Table 2. Conditional fidelities of the 5 individual temporal modes, each encoded as 
|D). Shown are also (for earch mode) the average number of photons g, in storage 
efficiency rj, the unconditional noise floor and the gi parameter (see main text). 


APPENDIX B: Results of the different encoded states 

We give here the measured conditional hdelities for differently encoded input states. 
These fidelities were used to calculate the process fidelity of the memory. We see here 
that the noise pn is not dependent on the input state. 


|'h*n) 

v{%) 

Pn (10-3) 

hi 

Fidelity (%) 

|H) 

3.3 ± 0.3 

9.3 ± 1.3 

0.28 ± 0.05 

84.1 ± 0.2 

|V) 

3.7 ± 0.3 

12.3 ± 1.5 

0.33 ± 0.05 

84.0 ± 0.1 

ID) 

3.6 ± 0.3 

10.1 ± 1.2 

0.28 ± 0.04 

85.5 ± 0.1 

|R) 

3.1 ± 0.2 

11.3 ± 1.6 

0.36 ± 0.06 

82.6 ± 0.1 


Table 3. Conditional fidelities measured for a mean input photon number /r=1.4±0.1 
prepared into various input states Shown are also the storage efficiency g, the 

unconditional noise floor and the gi parameter (see main text). 
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We here show the measured transmission coefficients for the five temporal input modes 
and their measured conditional fidelities. These numbers were used to calculate the 
criterion plotted in figure [^a) (black dashed curve). 


mode 

Transmission (%) 

Fidelity (%) 

1 

33.8 

97.2 

± 

0.4 

2 

28.0 

96.8 

± 

0.5 

3 

30.4 

97.4 

± 

0.4 

4 

30.1 

97.6 

± 

0.4 

5 

25.5 

97.0 

± 

0.5 


Table 4. Conditional fidelities of the individual transmitted temporal modes measured 
for an input encoded in |R) and g = 1.4. 


APPENDIX D: A quantum behaviour of the memory 

In this section, we derive the criterion that was used to infer a quantum behaviour of 
our memory, given that we probe it with coherent states only. 

We first summarize an already existing criterion [32], which only takes into account the 
non-unit efficiency of the memory in terms of the output. We then modify the criterion 
by also considering that a less than 100 % efficient memory cannot entirely absorb the 
input state. The main idea relies on the fact that the memory transmits part of the input 
state with a very good fidelity, so that a classical device that would like to reproduce 
the behaviour of the memory should also reproduce this transmitted state. 

4 . 1 . “Measure-and-prepare strategy” on the output only 

To infer a quantum behaviour of our memory, we propose to show that a particular 
protocol based on a measure-and-prepare strategy cannot mimic the experimental 
results. In other words, if we suppose that instead of the memory we place an 
eavesdropper, traditionally named Eve, that has only classical storage means at his 
disposal, we would like to know how well she could simulate the experimental results. 
It is known that with this strategy and an input qubit encoded into exactly n photons, 
the best achievable fidelity is 03: 

n -I- 2 

In our experiment we do not use Fock states with a well defined number of photons n but 
we use a coherent state with a mean photon number fi, so that equation (|^ has to be 
modified. In theory. Eve could take advantage of the coherence of the state to optimize 
her measurements on the state and maximize the fidelity of the states she re-sends. It 
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is though implicitly assumed that the experimentalist that prepares the state can easily 
blur the phase of the coherent state so that the stored state is only a statistical mixture 
of Fock states, encoding the polarization onto states of the form: 

jjTl 

Pin = y^ P{ti^,n)\n) {n\, with P{p,n) = e~^ — . (7) 

nl 

n=0 

Then, the output conditional hdelity is simply the statistical mixture of the hdelities 
for n ^ 1 photons: 


Fm{v) = 


+ 00 

E 


n=l 


n + 1 P{p, n) 
n + 2 1-P{p,0) 



e ^ — {1 + 

/i2 



• ( 8 ) 


Note that since we use photon counting techniques we do not measure the vacuum state 
|0), hence the sum in the equation above starts at n = 1. This is then by dehnition a 
conditional hdelity, i.e. conditioned on the detection of at least one photon. 


4 . 2 . Efficiency consideration 


In |32lll5], the authors have pushed the considerations one step further by taking into 
account the output inefficiency of the memory, so that Eve does not need to output a 
state each time. Eve can use this to only send an output state only if the number of 
photons is high enough in the input (n ^ Umin)- The conditional hdelity is then higher, 
and can be written as: 


Pm{p,Vm) 




+ 1 




+ 2 


7 


E 'PfPM 


n'^rirt 


n + 2 


7 




n^Urr 


( 9 ) 


where 0 < 7 < P{p, nmin) is a parameter that can be adjusted to mimic the memory 
efficiency rjM- 


4 . 3 . ‘Measure-and-prepare strategy” on the output and the transmitted states 
4 . 3 . 1 . Idea of the strategy 

In order to duplicate the complete behaviour of the memory, a classical scheme should 
reproduce its imperfections. One of these is the fact that the coherent state to be stored 
is not fully absorbed by the material, and part of it leaks into a transmitted state. We 
propose to implement a possible scheme based on the previously introduced classical 
measure-and-prepare strategies, shown on hgure 

In our analysis, we will suppose that whatever the strategy of Eve, she has to send 
a transmitted pulse which has exactly the same hdelity Ft and efficiency ry as the ones 
of the experiment, and the stored output pulse should be emitted with the same overall 
memory efficiency rjM SiS the one of the experiment. The parameter that she will try 
to maximize in the scheme is the hdelity of the output state Fm. The conclusion of 
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Figure 4. A classical scheme for imitating the quantum memory. The two measure- 
and-prepare strategies are shown in the two boxes and . 


a possible quantum behaviour of the memory will hnally be inferred by comparing the 
best achievable hdelity Fm using this strategy with the one obtained experimentally. 

To mimic the quantum memory with this classical scheme, we provide Eve two 
possible strategies: 

- She can, in a hrst strategy, measure the totality of the input pulse and try to infer 
the qubit state. She then sends her guess in the transmitted as well as the stored 
pulse. To correctly reproduce the transmitted and output efficiencies, a loss 6 can 
be adjusted in the output state. This strategy is depicted in the top part of hgure 

il 

- In a second strategy, instead of measuring the whole input state she can simply 
transmit a part of the state with a beamsplitter of transmission r] to reproduce the 
transmitted pulse, and measure the reflected state to reproduce the output pulse. 
As the hdelity of the transmitted state is generally not expected to be perfect, noise 
must be added to the state by sending a completely mixed state with probability 

q- 

These two strategies are complementary in the sense that the hrst one could lead 
to a too high output hdelity, at the cost of the transmitted state hdelity, while the 
second strategy could have the opposite ehect. Then, for each input pulse. Eve has 
the possibility to choose between these two strategies: the hrst one will be chosen with 
probability p, and the second one with probability 1 — p. 


4 . 3 . 2 . Mathematical derivations 
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Transmitted state: 

According to the previously depicted scheme, the fidelity and the efficiency of the 
transmitted state are: 


F, = - 


+ (1 - P)V 


Vt L 


l + q 


( 10 ) 

m = nu + 0 - ( 11 ) 

where Fm = Fm{p,Vm) and Fm is given by equation (N). As Ft and rjt are fixed in our 
scheme, the two previous equations can be inverted to find expressions for p and r/: 


p = 


V = 


Vt 


1±£ - Ft 
2 ^ * 


( 1 ) l+q _ p{l) 
'IM 2 n M 

( 1 ) 

Vt - PVm 


1 — p 


( 12 ) 

(13) 


Output state: 

Similarly, for the output state: 


Fm — — 
Vm l 


P^VmFm + (1 - p)(1 - v)VmFm 
Vm = pdVM + (1 -P)(l - v)Vm^ 

(O') fol 

with Ffj = Fm((1 — v)TiVm)- Given that rjM is fixed we have 


(14) 

(15) 

(l-p)(l-r/)' 

The strategy of Eve is to search for the maximal fidelity Fm by optimizing over the 
parameters d and q. 


( 2 ) 

vW 


Vm - pSpM 


Existence of a solution: 

We first show that a solution reproducing the measured efficiencies and the 
transmitted fidelity always exists in our model, which might not maximize the output 

(O') 

state fidelity however. If we set p = 0, g = 2Ft — 1, rj = rjt and rjf^ = tjm, fhe setup 


can fulfill all the conditions (12), (13) and (16), so that there is always a possibility to 


reproduce the conditional fidelity of the transmitted state Ft and the efficiencies rjt and 
rjM of the transmitted and output states respectively. It has to be noticed that in this 
case, only the second strategy is used (p = 0). 


Maximum ’’classical” fidelity based on our memory parameters: 

In figure we show the maximized conditional fidelity as a function of the mean 
photon number p of the coherent state, for Ft = 97.2%, rjt = 29.6% and rjM = 3.85%. 
The figure reveals that the maximum achievable fidelity with the classical strategy is 
lower by taking into account the transmitted state. 
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Figure 5. Maximum achievable conditional fidelity of the output state Fm by using 
the two measure-and-prepare strategies discussed here. Solid line: use of strategy 
presented in [32] (see eq. with rj = 3.85%. Dashed line: use of the strategy described 
here, for Ft = 97.2%, gt = 29.6% and rjM = 3.85%. 
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